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PROBLEM 


Develop the theory of internal waves in the ocean. Design and perform 
experiments to obtain needed military information. Specifically, investigate the 
diffusion of a solute through discrete layers of solvent, and develop a technique 
for measuring diffusion without disturbing the diffusion process. 


RESULT 


A technique was developed for determining the rate of diffusion between 
two discrete layers by taking minute samples at regular intervals with a microcell 
and determining concentration by conductivity measurement. 


RECOMMENDATION 


Apply the technique to investigation of the effect of wave action on dif- 
fusion between discrete layers. 


ADMINISTRATIVE INFORMATION 


Work was performed under SR 104 03 01, Task 0594 (NEL L41171) between 
January and September 1962. The report was approved for publication 26 July 1966. 
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INTRODUCTION 


Three approaches have been made to the study of internal waves and 
associated phenomena ~ theoretical analysis, model study, and field experiment. 
This report discusses the results of the first two applied to the specific problem 
of the diffusion of a solute between two discrete layers. The fluid model is pre- 
pared by dissolving different amounts of solute in successive layers of fluid. 


THEORY 


Attention is restricted to a two-layer liquid system with combined thickness 
]. Let the upper layer, of thickness h, initially be a pure solvent lying on a layer 
of thickness (J — h) of solution having a concentration Cy. If the diffusion coeffi- 
cient k can be taken as constant for concentrations ranging from 0 to C», the con- 
centration satisfies the simple diffusion equation, subject to the appropriate initial 
and boundary condition. This problem, as a special case of a more general treat- 
ment, is treated in the appendix. The solution is 


¢ fas sas 
e= SF fe( Zt) —ext( St _). ent( 22.) - ent ( Chi ) (1) 
2 VARt VAkt V4kt VAkt 
where erf stands for the error function, z is distance from surface, the absolute 
value of the argument is used, and the sign of the term is changed if the argument 


is negative. 


EXPERIMENTAL VALIDATION 


The experiment checks the results of the theory for a special case. If the 
agreement between theory and experiment is good, the theory and the assumptions 
on which it is based can then be applied with confidence to a range of similar 


problems. 
The basic problem is to get experimental curves. A sequence of measure- 


ments of concentration at several points along a vertical line is required, repeated 
several times, preferably at times varying exponentially. It is clear that the instru- 
ment used in making the measurement or extracting a sample should disturb the layer 
system as little as possible, both in moving into the fluid and in extracting a sample. 
The measurement of the chlorinity of very small samples of NaCl solution was con- 
sidered. However, as this method would have required many insertions of the sam- 
pling probe, and as adequate microchemical techniques were not available, it was 
decided instead to make conductivity measurements with a special microcell that 
allows a given sequence of measurements to be made with a single insertion. 

Since both the diffusion coefficient and the conductivity of the salt solution 
depend on temperature, a temperature bath was required. 


Apparatus 


MICROCELL AND BRIDGE 


The microcell was built into a capillary tube with an enlarged lenticular 
chamber. The microcell consists of two platinized platinum electrodes firmly 
spaced within a glass chamber with a volume of 0.01115 ml between its tip and a 
level % inch above the upper electrode. 

The bridge is basically an ac Wheatstone bridge excited with a low- 
voltage, single-stage vacuum-tube amplifier and equipped with a null-indicator 
tube. The bridge is provided with a vacuum-tube oscillator circuit supplying 1000- 
c/s bridge current. The operator balances the bridge by selecting a ratio and 
adjusting a resistance in the fourth arm of the bridge. A capacitor in parallel with 
this resistance can be varied to nullify the effect of the capacitance of the leads 
from the bridge to the cell. 

A Hamilton Model 0010 Micro-Syringe was used to withdraw the liquid 
through the conductivity cell. It consisted of a 1-ec syringe with a plunger in 
which a lead screw was mounted to withdraw minute amounts of liquid; as the lead 
screw was turned, the inner part of the plunger moved in or out, changing the vol- 
ume of the liquid in the syringe. The change in volume per turn of the lead screw 
was (0.0076 ml +1.3 percent. 

The lead screw could be withdrawn a total of 33 turns. In practice, it was 
turned three turns at each new position to insure that the sample was first thoroughly 
flushed from the chamber. Thus, 11 samples could be obtained for each insertion 
of the probe assembly. The assembly was securely attached to a feeding mecha- 
nism, a 12-cm, 28-threads-to-the-inch screw on which was attached a turn knob 
divided into 100 equal parts. The probe assembly could be advanced to each posi- 
tion with a least count of 0.01 mm. 

A cathetometer was used to determine depth (refer to Procedure). 


TANK AND BATH 


A 5-gallon tank contained the two-layer system. Temperature was main- 
tained at 21.00 (40.1) °C by immersion of the 5-gallon tank in a 15-gallon aquarium 
tank. The latter tank was insulated from the atmosphere by a 3-inch layer of 
Microlite, a commercial form of spun glass. A window cut into the side of the 
Microlite enabled measuring with the cathetometer, and a small aperture cut into 
the top admitted the microcell. Water maintained at 21.00 (40.02) °C flowed 
through copper tubing coiled around the 5-gallon tank. The temperature was regu- 
lated by a controlled-temperature bath manufactured by the Precision Scientific 
Company (fig. 1). 


Figure 1. Test setup, showing insulated aquarium tank, temperature control system, 
and conductivity cell. Photograph at right shows conductivity cell in entirety. 


Procedure 


The specific conductance bears the following relationship to the resistance: 
C=K,/R. The cell constant K, is determined with a solution of known specific 
conductance. In this instance, KCl with normality of 0.01000 was used at 21.0°C. 
The specific conductance is 0.001305 mho/em. The specific conductance of nine 
different NaCl solutions with normalities ranging from about 0.001 to 1 normality 
was measured; the temperature was maintained at 21.0 (40.1) °C. The cell constant 
did not actually have to be determined; the resistance of the cell corresponding to 
the normality would have sufficed. The results are plotted in figure 2. 

A sharp discontinuity in the concentration was achieved at the interface by 
laying a sheet of Dow Handi-Wrap over a layer of 1.05-normality NaCl. Distilled 
water at 21.0°C was poured over the sheet until the desired thickness was obtained, 
and the sheet was carefully withdrawn from between the two layers over a 5-minute 
period. Diffusion started at the time the trailing edge of the sheet passed the 
region in which the microcell would operate. The depth / of the two layers and 
the thickness A of the upper layer, together with the position of the microcell, were 
determined by means of the cathetometer. 

A short time after the layers were laid down, the first sequence of measure- 
ments of conductivity was made at points straddling the initial interface. The 
microcell was initially filled with distilled water and then inserted into the liquid 
to the first position, which was determined with the aid of the cathetometer. The 
screw knob on the micro-syringe was given three turns to withdraw 23 microliters 
of solution, on which the specific conductance was measured. The time was noted. 
The microcell was then advanced to the second position by means of the lead screw 
feed mechanism, and the measurement procedure was repeated. Eleven measure- 
ments were taken in a given sequence. The tip position was checked again with 
the cathetometer at the last position; then the microcell. was removed slowly from 
the body of the liquid, flushed out with distilled water, and made ready for the next run. 
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Figure 2. Calibration of conductivity cell. 
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RESULTS AND CONCLUSION 


Five sequences of readings were taken at approximately exponentially 
increasing times. For each sequence of readings, equation (1) in THEORY was 
used to determine a theoretical curve. The curve took into account the fact that 
an appreciable time was required to complete a sequence of readings. The nor- 
malized normality is plotted versus the normalized depth, with time a parameter 
along the curve. The conductance measurements at a given depth were entered 
into figure 2 to determine the normality, which was then plotted on the same figure 
as the theoretical curve (figs. 3-7). The agreement between theoretical predictions 
and experimental results is good, and it is concluded that the theory is valid, and 
can be applied to a wide range of similar problems. 


RUN 1 ¢ (AVERAGE TIME) = 23 MINUTES 


N, = 1.00 NORMAL 
k = 7.61 x 10-* CM?/MINUTE 


NORMALIZED NORMALITY (N/N,) 
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VERTICAL DISTRIBUTION OF CONCENTRATION 
(RATIO OF DISTANCE FROM SURFACE TO THICKNESS 
OF COMBINED LAYERS — 2/1) 


Figure 3. Run one -- diffusion profile at 23 minutes. 
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NORMALIZED NORMALITY (N/N,) 


RUN 2 ¢ (AVERAGE TIME) = 96 MINUTES 
N, = 1.00 NORMAL 
k = 7.61 x 10 CM’/MINUTE 
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Figure 4. Run two — diffusion profile at 96 minutes. 
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RUN 3 i (AVERAGE TIME) = 243 MINUTES 
N, = 1.00 NORMAL 
k =7.61 x 10-* CM?/MINUTE 
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Figure 5. Run three — diffusion profile at 243 minutes. 
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NORMALIZED NORMALITY (N/N,) 


RUN 4. ¢ (AVERAGE TIME) = 1330 MINUTES 
N, = 1.00 NORMAL 
k = 7.61 x 107* CM’/MINUTE 
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Figure 6. Run four — diffusion profile at 1330 minutes. 
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RUN 5 ¢ (AVERAGE TIME) = 3083 MINUTES 
N, = 1.00 NORMAL 
k = 7.61 x 107~* CM?/MINUTE 
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Figure 7. Run five — diffusion profile at 3083 minutes. 
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APPENDIX: ANALYSIS 


The analysis follows somewhat that in Sommerfeld.* 
The diffusion equation in one dimension is 


aus 1 du 


DE Os (1) 
Ozamivot 


where u is the concentration of solute, t is time, and k is the diffusion coefficient. 
Assume that the variables are separable; that is 


u = Zz) T(t) (2) 


2 
Use -(=) as the separation coefficient (where 2/ is the length of the specimen). 


Then: 
z+ (@) 2-07 + (“*) er =o (3) 
The solutions become 
TZ . [anz Sle kt 
La COS C=) b, sin ee) 5 Se (4) 


A general solution is obtained from the superposition of a number of such 
solutions for a proper set of values of n. 


Case of the Ring 


Take z = 0 at some point on the ring, and assume that the initial concen- 
tration f(z) is arbitrary but symmetrical about z = 0. The Fourier expression** 
for f(z) becomes 


< L(t 1 
fits) = y Anes (=), Ae = 3 i flz)dz, An = ae f(z) cos ("dz (5) 
n= 


The solution for u becomes 
Sc 
u(z,t) = » Ae \E cos(=) (6) 
n=0 


Now consider f(z) to be a ‘* function’”’ situated at z = +€. The expression 
for f becomes the theta function.*** Thus 


**The sine terms drop out because of symmetry. 
***Tt is usually given for € = 0 and for / = 1. 
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0(z,€;t) + O(z,-€;t) 1 . ~(42)'ne (z-€) 
u(z,t) = ep pal Se) | 5) ee, » a - cost ES) 
21 21 l 
n=| (7) 
= (nt) pe 
ate heed (7) eee "| 
n=] 


Application to a Tank 


Here instead of the ring, consider a rod of length / with the boundary 
conditions 


OB 0; z=0,1 (8) 

Oz 
expressing the fact that there is zero flow of solute at the ends. If the rod is 
reflected at z = 0, with the reflected part having the same initial distribution of 
solute and the free ends joined, it is seen that the boundary conditions (8) are 
automatically satisfied, and the problem is reduced to that of the ring. Thus, 
immediately 


1 f(€) 


or 
i ie) oy ae eos + ena cay (AS d&~ (10) 
7Jo . i i v 
ite 


The theta function given by expression (7) does not converge rapidly for 
very small values of kt. We shall transform it to a very rapidly converging series, 
but first we must put it into a suitable form. Consider 


1 = -() mel 
8(z,+€ 5t)) = Ak mY) » oN Cegeen a any 
n=] 
Let 
jo? Xs > oy (12) 
Then 
CMDES FU) Zila 9 Sy da" cos Qnnx, | = os 8(y ,|7) (13) 
r au 1 an 21 a 
n= 


The famous transformation formula for the theta function gives 


-¥% =—wiy? 
O(y|T) = () IQ U6 o(% -*) (14) 
1 


Substituting (12) into (14), (14) into (13), and (13) into (10) gives 


l on l(2-§) +2nl}?/4ke , po l(2t&) + 2nlP/4ke 
(z,t) = d il 
u(z,t) = Ih f(§) 7 Se (ae Es) 


—oo 


Each term represents a concentration pole at €+2nl or —€+2nl, n= ...., —2, 
—1, 0, 1, 2, ...., and thus the initial concentration at € leads to an infinite sequence 
of concentration poles made up of the primary pole and all its images reflected at 
z= 0,1. The primary heat pole 


—(z-§&)? 
Ln 
4 kt 


could have been deduced initially (it is called the principal solution), and the 
images could then have been introduced to satisfy the boundary conditions. 


Special Case 


Suppose that initially in the tank there is a distribution of concentration 


f(§) = 0, O0<z<h 


(16) 
i) so, WK<a<t 
Then (15) becomes 
== —_|/! > eo ey on (2t§+2nl)/4kt| ge (17) 
ae t 
By letting 
(z+ §+2nl)? : 
Sete Mie EE Ges 1 
Fe r (18) 
we transform it into 
a io) z—hta2nl ztlt2nl 
We) a2 = Vise omo=la's 2 Vie Wem ds (19) 
2 Na z—1+2nl < Vir thee 
ae V4kt 
The error function is defined 
2 (xX 2 
ere Cmmids (20) 
Tear iF 


Ug) 


20 


Note that 
2 : ; 
=|" oP tg = |* BO? Fe = f° OO Gls aor) ora) 
Vr’ a Viry—° Vn 


The error function is defined for a positive upper limit of integration. Here the 
upper limit may be negative. Note that 


2 oe ds = 1—erfa, a>0 
Ji I~ 
Thus 


ulz, | Dy, ie ee ae : IE a) BEf Cafes 


where the sign of the term is changed if the argument is negative. 


As an example of the special case, consider a layer of solvent of thickness 
h such that h/] = 0.2. Let this layer initially lie on a layer of solution with con- 
centration C,. The subsequent distribution of concentration C divided by C, 
(that is, C/C,) is plotted as a function of z// for the following values of the non- 
dimensional parameter kt//?: 0.01, 0.0004, and 0.0001 (fig. Al). 


The time required to reach the above indicated distributions from the initial 
state is given in table Al. Here / = 50 cm, solutions are aqueous, concentration is 
in gm-molecules per liter, and k is given in cm’/day. 


TABLE Al. EXAMPLES OF DIFFUSION TIME 


ve Diffusion 
Coefficient (k) 


Material | Concentration | Temperature, Time (t), days 


It can be shown that the gradient at h// = 0.2 is (for © not too large) 


proportional to 1/Vkt and directly proportional to J. For a given gradient, the time 
required for a sharp interface to deteriorate to that gradient is inversely proportional 
to k. 


kt/l? = 0.0001 


kt/I? = 0 


C/C, 


Figure Al. 
of time. 
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Vertical distribution of concentration as a function 
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